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Abstract 

The Social Force Model is one of the most prominent models of pedestrian dynamics. As such naturally 
much discussion and criticism has spawned around it, some of which concerns the existence of oscillations in the 
movement of pedestrians. This contribution is investigating under which circumstances, parameter choices, and 
model variants oscillations do occur and how this can be prevented. It is shown that oscillations can be excluded 
if the model parameters fulfill certain relations. The fact that with some parameter choices oscillations occur 
and with some not is exploited to verify a specific computer implementation of the model. 


1 Introduction 

The Social Force Model of pedestrian dynamics is a model that aims at describing the movement of pedestrians 
with the predominant purpose of simulating pedestrian movement on computers. The force of pedestrian /3 on 
pedestrian a typically has the form 

fal3 = eaf} ( 1 ) 

where ^() is a function which grows with increasing distance between both pedestrians and can depend on the 
velocities of one or both pedestrians. The function ie() suppresses forces the more pedestrian 13 is located outside 
the current walking direction of pedestrian a. 

The Social Force Model has first been introduced in 1995 [1]. This variant later was called “elliptical speci¬ 
fication I”. A second variant (circular specification) has been proposed in 2000 [2] and a third variant (elliptical 
specification II) in 2007 [3]. The difference between the three variants lies mainly in the way the velocities of 
two interacting pedestrians are considered in the computation of the force between them. The 1995 variant 
considers only the velocity of the pedestrian who exerts the force. The 2000 variant does not at all consider 
velocities (only the distance between pedestrians) and the 2007 variant considers the relative velocity between 
both pedestrians (the pedestrian who exerts the force and the pedestrian on whom the force acts). For the 
analytical considerations in this paper mainly the simplest variant from 2000 will be considered. Nevertheless it 
will also be discussed how results will change qualitatively if the variants as of 1995 or 2007 are applied. 

Under “oscillations” in this paper unrealistic artifacts in the trajectory of a pedestrian approaching another 
pedestrian, his destination or a wall is understood. The occurrence of oscillations in the sense of this contribution 
has been discussed in a number of contributions [4-7] and it is often claimed that oscillations cannot be avoided 
in the Social Force Model, but that they just can be made small. In this paper it will be shown that this is not 
correct and exact conditions for the value of the model parameter such that oscillations occur will be derived. 

In the remainder of the paper first a single pedestrian approaching a destination is investigated, then a 
pedestrian approaching another pedestrian who is standing still and finally two pedestrians approaching each 
other. In each case the model is reduced to one dimension and the noise term is set to zero. In the first section 
on approaching a destination the problem will be shown to be most severe as with certain conditions oscillations 
cannot be prevented and continue infinitely long. At the same time - as will be argued - for this case it is not 
very relevant, as there are simple, pragmatic solutions. The second and third case yield restrictions to the choice 
of parameters which can produce realistic, oscillation-free behavior. 


2 A pedestrian approaching a destination 

In this section we are interested in and discuss the equations of motion and their solution of a single pedestrian 
approaching a destination coordinate (i.e. a point) where he is required to come to a standstill. We assume that 
in the beginning the pedestrian is walking with his desired speed vq straight towards the destination coordinate, 
so there is no tangential component of the walking velocity. Then we can describe the pedestrian as walking 
from positive x coordinate into negative x direction towards the destination which is at x = 0. Since the 1995, 
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2000, and 2007 variants of the Social Force Model only differ in the force between pedestrians and not the driving 
force term all results of this section hold for all three variants. 

We assume for now, that the desired velocity is always some externally given vq and is always pointing from 
the pedestrians current position towards x = 0. This assumption is the simplest one and it can be questioned - 
as we will do below. With it it is obvious that there will be oscillations around x = 0. Our intention here is to 
investigate the quantitative details of these oscillations. 

In general the equation of motion for this pedestrian reads 


x{t) = 


—sign{x{t))v{) — x(t) 


( 2 ) 


where r is an external parameter which typically has values between 0.1 and 1.0 seconds. 

We require the pedestrian not only to reach x = 0, but also to stand still there as arrival condition. Because 
the pedestrian has a speed larger 0 (or, considering walking direction: smaller 0) he will walk over the destination 
and be on the left (negative) side of x coordinates. There the desired velocity points into the direction of positive 
X coordinates. So we have for the time following the moment when the pedestrian is at x = 0: 


x{t) = 


■yo - x{t) 


(3) 


This is solved by 


x{t) = vq — ae 

_^ 

x{t) — 6 + xot + are ^ 


(4) 

(5) 


where a and h are integration constants which need to be determined by initial conditions. 

We choose t = 0 at the moment when the pedestrian is at x = 0. Then x{t = 0) = —vq. However, for later 
usage we want to set here more general x(t = 0) = —u and remember that for our particular case u = vq. With 
the two conditions x(0) = 0 and x(0) = —u we can determine the values of the integration constants: 


a = Vo u 

b = — (xo + i^)t 


So we have 


x(t) = Vq - {vo + u)e ^ 
x{t) — vot — (xo + u)t ^1 — e~ 


(6) 

(7) 

(8) 

(9) 


Now we can compute the time ttum when the pedestrian stops (and turns around) x(to) = 0 and the position 
x(to) at which this happens: 


x(t 


tturn 

1 Xo + X 

= T In- 

Xo 

'turn) 

= TXo ^In ^1 + ^ 


u 

Vq 


( 10 ) 

( 11 ) 


In the initial case, when u — vq this simplifies to Uum = Tln(2) and x{tturn) — TXo(ln(2) — 1). 

This is only half the way to go. The actual question is how fast a pedestrian returns to x = 0 when he has 
passed it before with speed u and how long such a loop takes. 

Now we choose t = 0 for the moment when the pedestrian is standing still at x{tturn)- The time tretumed 
at which the pedestrian returns to x = 0 therefore has to be understood as time after Uum and not as absolute 
point in time. 

We begin again by determining the value of the integration constants. From equation (4) and x(0) = 0 
follows a = xq. With equations (5) and (11) we have 


x(0) = b-\- vqt = xoT ( In ( 1 + — ) - — 

' ' Vq J Vq 


b = XOT ( In ( 1 + — ) - ( 1 + — 

Vq J V ^0 

.(0 = -or(ln(l + M-(l+j) + t+e- 


Because x(tretumed) = 0 wc havc the following equation for tretumed- 

n ^returned , - ^returned 

0 = xor ln lH-- IH-H-he 

V V '^^0 / V ^0 / T 

0 = (preturned “ « + ln(a) + 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
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with the substitutions (j^retumed = ^returned It Uud a — 1 + u/ Vq. Wg do another substitution Lp — 4>returned — 

a + ln(a), transforming the last equation to 


0 = ^ + (17) 

—ae~^ — (18) 

Obviously one solution is p — —a, but the general solution is 

p^W {-ae-^) (19) 

where W{) is the Lambert W function [8-10], which by definition is the inverse relation of f{y) = ye^. 
Resubstituting p and (j^retumed (for reasons of convenience we do not resubstitute a) we have: 

^r5^^=lK(-ae-“)+a-ln(a) (20) 

In the interval y G [—l/e..O] W(y) has two branches, denoted W-i{y) and Wo(2/); and with 1 < a < 2 for 
—ae~^ we are within this interval where W{y) has two branches. It holds 

lL_i (-ae“^) = -a (21) 


In this case it would be tretumed = —ttum- Thus the W-i branch gives the backwards in time solution which 
we are not interested in here (because we already have computed it above). Therefore we have to continue with 
the Wo branch for which 

ITo / -a (22) 

although W is the inverse relation of f{y) = ye'^. In the remainder we write W for Wq. 

In case u = vq, i.e. a = 2 the numerical value of the solution is + 2 — ln(2) = 0.90047. 

Using equation (20) in equation (4) we get the speed x{treturned) at the time when the pedestrian returns to 
(r = 0 in dependence from which is the speed at the time when the pedestrian last was at x = 0: 

X{treturned) = Uq (l + IT (23) 

where we have used the defining equation of the W function y = W{y)e^^'^\ 

From here on the properties have an index that determines the recurrence to the offspring. For example to 
is the (absolute) time when the pedestrian is for the first time at the offspring, ti denotes the (absolute) time 
when the pedestrian returns for the first time to the offspring and so on. In the case of properties that do not 
describe passage of the offspring, but the loop durance and turning position the index enumerates the loops. 
This means that for these two properties there is no index 0 and that = tn-i + Atn- 

Equation (23) means that for the n + 1th passage of the offspring an+i depends on an like 

an+i = 2 + W {-ane-’^^) (24) 

The time At^+i it takes for the n + 1th loop depends on an like 

=a„ + H^(-ane-“"). (25) 

Rewriting equation (11) in terms of a and writing \xn\ for the turn around distance of the nth loop we have 

= ttn - 1 - In(an) (26) 

TVo 

Table 1 shows the results from the first 30 passages if the process begins at t = 0 and with speed vq. Figures 
(1) to (3) visualize these data. 

As typical values are r = 0.4 s and vq = 1.5 m/s it takes about 7 passages before the amplitude \xn\ gets 
smaller than 1 cm. This takes about 2.1 seconds. The 7th oscillation has a period of 0.15 seconds. To resolve 
this a computer implementation would have to have a time step which is at the very maximum half that value 
0.075 seconds or 14 simulation steps per second. 

Because an+i = 1 only if W{—ane~°''^) = —1 which in turn is only the case if an = 1 there will be infinitely 
many oscillations and in A a proof is given that the sum of the Atn diverges, i.e. the pedestrian will oscillate 
infinitely long around the destination coordinate. 

At first sight this is an unsatisfactory result with regard to the Social Force Model, as such oscillations are 
unrealistic. However, we have to ask how realistic our assumptions and initial conditions were. This concerns in 
particular the desired velocity vq. We demanded in the beginning that the pedestrian should come to a stop at 
X = 0. Nevertheless we have set the desired velocity all of the time to one particular value |uo| > 0. This is too 
simplistic. Real persons plan ahead and adjust their desired speed: if they desire to stop at a certain position 
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Passage/Loop(n) 

l^n 1 
TVn 

T 

tn 

T 

Un 

vn 


0 


- 

0.000 

1.000 

2.000 

1 

0.307 

1.594 

1.594 

0.594 

1.594 

2 

0.128 

1.018 

2.611 

0.424 

1.424 

3 

0.071 

0.754 

3.365 

0.330 

1.330 

4 

0.045 

0.600 

3.966 

0.270 

1.270 

5 

0.031 

0.499 

4.465 

0.229 

1.229 

6 

0.023 

0.428 

4.893 

0.199 

1.199 

7 

0.017 

0.374 

5.266 

0.175 

1.175 

8 

0.014 

0.332 

5.599 

0.157 

1.157 

9 

0.011 

0.299 

5.898 

0.142 

1.142 

14 

0.005 

0.199 

7.062 

0.096 

1.096 

19 

0.003 

0.150 

7.898 

0.073 

1.073 

29 

0.001 

0.100 

9.088 

0.049 

1.049 


Table 1: Numerical results for the first 30 passages resp. oscillations. The value for is understood before an 
oscillation, while is the period of that oscillation, and the total time t^jr after that same oscillation. All 

values are dimensionless. The numerical values of the Lambert W function have been computed using Wolfram 
Alpha [11]. 


oscillation no. 


Figure 1: Left: Value of an at the beginning of oscillation n. Right: Maximum amplitude \xn\/{vor) of oscillation 

n. 






oscillation no. oscillation no. 

Figure 2: Left: Period Atnir of oscillation n. Right: Total time tn/r at nth passage of offspring (i.e. time after 
nth oscillation; t=0 is when the pedestrian reaches the offspring for the first time and with a speed vq). 
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Figure 3: This shows the evolution of the amplitude Xn/{voT) over time in [s]. 


they adapt their desired velocity beforehand to just achieve that. So we have to ask how we have to modify vq 
dynamically to account for that. 

More precisely we ask: at which distance db does a pedestrian have to start braking in the Social Force 
Model, if he sets Vb as desired speed opposing his current speed u7 And how long does it take before he comes 
to a stand still? 


For this the equation of motion reads 

x{t) = 

-x{t) 

T 

(27) 

with the following (initial) conditions: 

x{t = 0) 

= —u 

(28) 


x{t = 0) 

— d 

(29) 


x{to) 

= 0 

(30) 


x{to) 

= 0 

(31) 


where we have defined the moment where the pedestrian starts to brake as t = 0 and the moment where he stops 
at X = 0 as to; d is the distance before stand still at which braking has to begin which we are looking for. 
Solving this results in 


to = rln^lA^^ (32) 

d, = r.,(j-ln(l+J)) (33) 

This gives finite and positive to and db for all r'b > 0. Thus, it is not sufficient to set the desired speed to zero 
for a pedestrian who wants to stop as this would imply an infinitely large braking distance and an infinitely long 
time to come to a stand still. If we assume that a pedestrian for braking can at maximum have a desired speed 
Vb = Vo the minimum time for braking is to = Tln(2) and the minimum braking distance is db = rvo{l — ln(2)) 
to come to standstill from an initial speed u = vq. 

A different and pragmatic solution of the problem is that in real-world planning applications usually there is 
no requirement set for the speed at which a pedestrian arrives at a destination. It is furthermore usual that as 
destination not a point is given but that a cross section has to be crossed or the pedestrian has to move into a 
specific area. If no restriction is given for the speed at arrival the problem disappears. 

A third objection is that also a real person cannot comply to the request to stand still with his or her center 
exactly at some given coordinate. Real people always sway a little bit when they (try to) stand still [12]. Why 
then should one require this from a simulated pedestrian? The oscillations around the destination point which 
were found here in their functional form surely do not match the swaying of real people when they stand “still”, 
but the amplitude of the ones of the model quickly fall below those of real people. When the real swaying is not 
required to be reproduced by models implicitly a required limit to precision is set and any movement of that 
order of magnitude or below should be acceptable. 
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3 A pedestrian approaching a standing pedestrian 


3.1 Theory 

Now we want to investigate the situation when some pedestrian M is approaching another (pedestrian S) one 
who is standing at a fixed position - again in one dimension, for an empirical investigation see for example [13]. 
Since the pedestrian who is exerting the force is not moving, the 1995 variant (elliptical specification I) of the 
Social Force Model will produce the same result as the 2000 variant (circular specification) which is the one we 
are investigating, but different from the 2007 variant (elliptical specification II) where the speed of the pedestrian 
on whom the force acts modifies the strength of the force. 

Assume pedestrian S is standing still at a: = 0, facing into positive x direction and pedestrian M is approaching 
from there, i.e. pedestrian M has positive x coordinate and a speed directed towards negative x^ i.e. x(t = 0) < 0 
as well as the desired speed i;o > 0 which is constant in time. It is assumed that pedestrian M at time t=0 is far 
away from pedestrian S and walking with his desired walking speed x{t = 0) = —vq. 

Then the equation of motion for pedestrian M in the 1995 and 2000 variants of the Social Force Model is 


x{t) = 


—Vo — x{t) 
T 


X ( t ) — 2R 

+ Ae B 


(34) 


with R as radius of a pedestrian (we assume here for simplicity that all pedestrians have the same radius) and 
T, A, and B parameters of the Social Force Model. 

Pedestrian M will come to a standstill at distance 


ds = Bln 



+ 2i? 


(35) 


As pedestrians should not bump one into another ds must be larger than 2R. Therefore we have a first condition 
for the parameters of the Social Force Model to yield realistic results: 


At > Vo 


(36) 


We assume that oscillations only occur generally if they also occur close to equilibrium. Thus, if there are 
no oscillations when pedestrian M is already close to x = ds then there are no oscillations at all. Expanding the 
exponential function in equation (34) into a Taylor series around x — ds gives 

m + ^ + « 0 ( 37 ) 

Equation (37) is the equation of the damped harmonic oscillator, with the known three solutions: under damped 
(pedestrian M approaches the equilibrium point oscillating around it), critically damped (pedestrian M ap¬ 
proaches the equilibrium point as quick as possible but without oscillations) and over damped (pedestrian M 
approaches the equilibrium point slower than would be possible without having oscillations. Which of the three 
cases is realized depends on the relation of parameters: 


B 

> 1 ^ under damped 

(38) 

4^ 

B 

= 1 ^ critically damped 

(39) 

B 

< I o over damped 

(40) 

Thus in addition to equation (36) with 

4:Vor < B 

(41) 


we have a second condition for the parameters of the 1995 and 2000 variants of the Social Eorce Model to have 
a chance to yield realistic results. 


3.2 Implications for parameters found in literature 

In the literature one can find values for some or all of the parameters A, B, r, and vq gained from empirical 
observations or laboratory experiments. If all four parameters are given, one can use relations (36) and (41) to 
test if realistic results for the scenario discussed here can be expected. If not for all four parameters values are 
given those two equations allow to set limits. Table 2 shows a compilation of such data and how they relate 
to or with relations (36) and (41). It can be seen that where all four parameters were taken from experiment 
relation (36) easily holds while relation (41) is clearly violated. Whereas if parameters A and B (and r) have 
been calibrated or given in a work then relations (36) and (41) require pedestrians to walk quite slowly, while if 
Vo and r are calibrated or given it needs relatively large values for A and B that relations (36) and (41) hold. 
This is an indication that with the circular specification (alone) the parameter space that yields realistic results 
is small. 


6 




Source 

A in [m/s‘^] 

B in [m] 

r in [s] 

uo in [m/s] 

eq. (36) 

eq. (41) 

[3] 

0.42 ± 0.26 

1.65 ± 1.01 

T > Uq /Arfjiax 

< 0.67 






^ ^max ij^ 




[2] 

26 . 67 t 

0.08 

0.5 

0.8 

OK 

violated 

[14] 

> 2.25 

> 3.34 

0.61 

1.37 



[15] 

0.16 ± 0.01 

4.16 ± 0.65 

0.5 ^ 

< 0.087 



[15] 

0.45 ± 0.41 

13.5 ± 18.4 

0.5 ^ 

< 0.43 



[16]' 

12.0 ± 0.2 

0.16 ± 0.08 

1.09 ± 0.35 

1.34 ± 0.21 

OK 

violated 


Table 2: Bold face marks value from literature; normal face is computed with equations (36) and (41). Where one 
parameter was calculated from another one where a range is given the range was utilized such that the parameter 
space for the derived parameter is as large as possible. from 2000N/75kg. Assumption and input for 
calibration. The simplified Nomad model discussed in said publication is identical to the circular specification 
of the Social Force Model, except that the radius of pedestrians is not stated explicitly in the exponent. If - what 
is not stated explicitly in said contribution - “distance” for the Nomad model means body surface to body surface 
the parameter meanings are identical, while if it means center to center distance the value of parameter A would 
change by a factor and even equation (36) could be violated. 

3.3 Elliptical specification II in comparison 

For elliptical specification II a rigid analysis is much more difficult than for the circular specification since the 
elliptical specification II does not only depend on the distance of the two pedestrians but also on their relative 
velocity. Still one can estimate the effect of the added consideration of relative velocity on the occurrence 
of oscillations. In the elliptical specification II of the Social Force Model [3] the force from pedestrian [3 on 
pedestrian a is defined as 


fcx{rcx,rp,rcx,rp) 


oc^idoc^ 1 'f'ci)^Q{da.j31 doiP^ 


(42) 



A 1 /I \ ^ ^ COs(^q:^ (da^; 

“T W ) 2 

.fa)) 

(43) 

cos(6>a^(dc/3,rk)) 


r OL ' docj^ 

\'^ol \ 1 


(44) 

doil3 


Fa - F/3 


(45) 

di^doc ^, d(xp^ 


Fq,/ 3 (do;^ 5 da^) 


(46) 

^oiisi^dcip 1 d(xp^ 


^a/3 Fq/3) 

AcxBoiG- 


(47) 

boc^i^doi^^ d(xp^ 


2 (|da^ 1 T |da^ T doL^Atoi |)^ 

- (dapAta)'^ 

(48) 

where r gives the position of pedestrians and A, B, A, and At are model parameters. Pedestrians’ positions and 
derived properties are time dependent, other properties are constant. In one dimension and with pedestrians 
facing each other this simplifies to 

COs{0cl3) 


1 


(49) 

W{0ap) 


1 


(50) 

doc^ 


X(x — XjB w.l.o.g. assumed to be > 0 

(51) 

foi {X(x 5 X(3 •) Xoif ) 


9 ida.P 1 da.pi) — , , ^apidap 

ddap 

1 dap) 

(52) 

^CX^{dcx^ ■) dcK^) 


^oil3(^a/3 ^^oip) 

AaBaG 


(53) 

bci/3{dcx/3 ^ daf3^ 


2 ^(da/3 4~ |da^ T daP Ata\)‘^ 

- {dapAta)"^ 

(54) 



0 for [dap + dapAta) < 0 


(55) 



\J^ap + dapdapAta othcrwisc 

(56) 
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( 57 ) 


Therefore the force is either zero if ((ia /3 

Xp, Xoi-i 


+ dcisAta) < 0 or it is 

— 9{d(y.f3 ■) dcxl3^ 


_ ^ ‘^doi^ T doc^^toi 

— ^ - 6 

^ Y docjsdcxjB ^ioL 

_ . da -i 

— Aoc^e 


^Ctl3~^^cx.^^cxl3 


(58) 

(59) 


with da being the arithmetic and dg the geometric mean of current and projected distance (da /3 resp. da^ + 
dai 3 ^ta). It can directly be seen that for large distances equation (58) reduces approximately to the circular 
specification from [2] which depends only on distance. 

For a pedestrian a approaching another pedestrian fd from positive x values d is positive and d is negative. 
Therefore and because of the inequality of arithmetic and geometric means it holds d^p + da/s^ta < dg < da < 
da/ 3 - Then obviously as long as d > —dAt the force as in equation (58) is larger compared to the circular 
specification and consequently pedestrian a will not overshoot for certain parameter choices where he does so 
with the circular specification. 

In case pedestrian a overshoots over the equilibrium point and turns around it would be desirable that the 
force is not larger but smaller than with the circular specification. However, since in this case d becomes positive 
it holds that da /3 < dg < da < da /3 + da /3 Ata and therefore in equation (59) the exponential factor gives a smaller 
value than with the circular specification, yet the fraction factor has a value > 1 and may for large values of 
parameter B outweigh the damping effect from the modification in the exponential function. There are three 
indications that also in this phase of pedestrian a’s movement and in general oscillations are suppressed with 
elliptical specification II: first, for large values of parameter B already in the circular specification there are no 
oscillations as equation (41) tells us. This means that where in an isolated view on the “way back” the problem 
is most pronounced the system may actually not even evolve to the point that it exists. 

Second, one can expand equation (58) in a series for small values of da^Ata/da^: 


f(xiXoii X ^ Xoi^ X ~ AcxC 



da.13 dfxj^At a. 
Bex dot 13 


(60) 


In the moment pedestrian a turns around it is da^ = 0 and therefore circular and elliptical specification II 
yield an identical force. Starting to move backward with now positive da /3 equation (60) tells us that elliptical 
specification II yields smaller forces with the difference to circular specification leveled for B^ ^ oo. 

Third, with the - admittedly arguable - additional assumptions that 

dap Atex dexp (bl) 

dexfsAta « Boc (62) 

one cannot only expand for small ^(t), but also reduce the complexity of equation (58) with regard to da/ 3 . 
Various approximate forms of that equation can be derived in this way of which one is analytically solvable and 
contains parameter Ata (omitting the indices a): 

i{t) + i (l + ^) C(i) + = 0 (63) 


where in comparison to equation (37) just the factor before ^(t) has changed. This leads to a less strict require¬ 
ment for avoiding oscillations: 

4:VoT < H + 

than equation (41). 

Finally a note on the case d + dAt < 0 where the force is zero: if pedestrian a starts to approach /3 from a 
sufficiently large distance and with a typical pedestrian walking speed it will be d + dAt > 0. Since equation 
(58) diverges to positive values at d + dAt 0^ the pedestrian will slow down and eventually be at rest or turn 
around before d < —dAt. For a departing pedestrian a it is d > 0 and thus always d + dAt > 0. Only when a 
simulation is initiated with d + dAt < 0 it may yield dynamics that do not fit into this line of argumentation; 
compare [17]. Such extrinsically prepared “Garden of Eden” states can be dealt with for example by a model 
extension (if one does not simply want to exclude them by construction). 


voAt\‘^ 

^ J 


(64) 


3.4 Simulations 

For realistic applications it makes sense to choose the parameters such that no oscillations occur. However, to 
verify a computer implementation of the Social Force Model it can be interesting to use parameters just around 
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the critically damped conditions and check if oscillations do and do not occur according to expectations. Model 
specific validation work can amend validation tests like those in the RiMEA test cases [18] which usually are 
formulated model independently. In this subsection such a model specific verification process will be carried out 
exemplarily utilizing PTV Viswalk [19]. 



Figure 4: Simulation scenario. 

Figure 4 shows how the one-dimensional setting has been implemented in the two-dimensional simulation 
model. The walking area is modeled to be 0.55 m wide. A pedestrian facing into positive x direction is stopped 
by a red traffic signal at x=0. Then a second pedestrian is set into the model at x=52 m moving into negative 
X direction. 


r 

distance 

expected 

difference 

0.7 

0.4556 

0.4570 

0.0014 

0.8 

0.4821 

0.4837 

0.0016 

0.9 

0.5044 

0.5072 

0.0028 

1.0 

0.5269 

0.5283 

0.0014 

1.2 

0.5645 

0.5648 

0.0002 

1.5 

0.6073 

0.6094 

0.0021 

2.0 

0.6655 

0.6669 

0.0014 

3.0 

0.7483 

0.7480 

-0.0003 

4.0 

0.8044 

0.8056 

0.0011 

5.0 

0.8482 

0.8502 

0.0020 


Table 3: Stand still distance [m] in dependence of parameter r [s], with A = 1.6 mjs^^ B = 0.2 m and vq = 1.5 
m/s. Column “expected” shows the values according to equation (35). 

Parameter settings: since A_socjnnean controls the contribution of elliptical II specification it has been set 
to zero (in the remainder ”A refers to A_socjiso in the software; for the pedestrian on the left it has been set to 
zero as well); A = 1.0 to minimize effects from the in fact two-dimensional nature of the simulation; stochastic 
noise has been set to zero . The value of the noise parameter has been set to zero for both pedestrians to not 
unnecessarily blur the results. 


B 

distance 

expected 

difference 

B 

dist. 

exp. 

diff. 

0.1 

0.5831 

0.5847 

0.0016 

4.0 

3.2870 

3.2880 

0.0010 

0.2 

0.6524 

0.6540 

0.0016 

6.0 

4.6734 

4.6743 

0.0008 

0.3 

0.7240 

0.7233 

-0.0007 

9.0 

6.7530 

6.7537 

0.0007 

0.5 

0.8590 

0.8620 

0.0029 

12.0 

8.8326 

8.8332 

0.0005 

1.0 

1.2073 

1.2085 

0.0013 

18.0 

12.9917 

12.9920 

0.0003 

2.0 

1.9005 

1.9017 

0.0012 

24.0 

17.1509 

17.1509 

0.0000 


Table 4: Stand still distance [m] in dependence of parameter B [m], with A = 2.0m/s^, r = 1.5 s and vq = 1.5 m/s. 
Column “expected” shows the values according to equation (35). The difference between theoretical expectation 
and simulation is in all cases below 3 mm. 

At first we investigate where the second pedestrian comes to rest. According to equation (35) this depends 
on the values of the parameters B_socJso, A_socJso, vq, t, and R. The latter was set to be 0.2577 m. Keeping 
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A^socjIso and ^'o constant and increasing the value of r in steps of 0.1 s the second pedestrian does not pass 
through the first one for the first time at a value r = 0.6 s and both do not overlap visually for r > 0.8 s. At 
T = 0.9375 s when At = vq the distance of the central points of both pedestrians is 0.5135 m which comes 2 
mm close to 2R. Table 3 and figure 5 show values for the stand still distance in dependence of some values for 
parameter r with all other parameters kept constant. The theoretical expectation is met well in all cases. 

Table 4 and figure 5 show values for the stand still distance in dependence of some values for parameter 
B^soc-iso with all other parameters kept constant. The theoretical expectation is met well in all cases. 

The stand still distances of table 4 are unrealistically high for all but the smallest values for parameter B. 
We have chosen the parameters in this way not to scan for realistic parameter values, but because with these 
parameters one can well demonstrate that for certain cases (small values for B) there are oscillations and in 
others (large values for B) there are none. Figures 6 and 7 show the time evolution of the position of the 
approaching pedestrian for various values of parameter B. 




Figure 5: Left: Visualization of the data of table 3. The expectation for the regression curve would be ^ = 
0.2000 ln(x) + 0.5283. Right: Visualization of the data of table 4. The expectation for the regression curve would 
be ^ = 0.6931X +0.5154. 




Figure 6: Left: Position of approaching pedestrian over time for various values of parameter B. Right: As A = 2.0 
mj^ vq = 1.5 m/s, r = 1.5 s the system is critically damped for 5 = 9.0 m. For increasing B the oscillations get 
smaller and vanish for B = 9.0. 


Neglecting for the under damped case the damping exponential function the approximately expectated time 
distance Tr between two reversal points in the under damped cases is: 


Tr = 



(65) 


Table 5 shows a comparison between actual and expected Tr for various values for B which generate under 
damped behavior. 


3.5 Two pedestrians approaching each other 

If two pedestrians approach each other in one dimension they may come to a stand still at an equilibrium distance 
or they might jointly (and with an equilibrium distance) move into one of the two possible directions. If all 
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Figure 7: Left: Position of approaching pedestrian vs. time for under damped and critical cases with regard to the 
value of B. Right: Zoom to the region of oscillations. 


parameters (vq, A, B, and r) are identical for both pedestrians one can rewrite the coupled equations of motion 
of both pedestrians as one equation for the movement of the center of mass and one equation for the relative 
motion. Latter one can again be approximated by an oscillator equation. Compared to the case above one has 
an additional factor 2 at the friction term implying that the system is critically or over damped for 

Svqt < B ( 66 ) 

which leads to the conclusion that the case of two mutually approaching pedestrians sets a stronger restriction 
on parameter choice than the case when one pedestrian is approaching another who is standing, at least if 
oscillations are to be avoided. The 2007 variant of the Social Force Model in this case suppresses oscillations 
even more than when one of the two pedestrians is standing still. 


B 

number of reverses 

simulation 

approximate expectation 

0.1 

11 

1.035 

0.999 

0.2 

7 

1.442 

1.421 

0.3 

3 

1.800 

1.750 

0.5 

3 

2.250 

2.286 

1.0 

2 

3.100 

3.332 

1.5 

2 

4.000 

4.215 

2.0 

2 

4.850 

5.038 

3.0 

2 

6.550 

6.664 


Table 5: Comparison of simulated and approximately expected time [s] between reversals. 


4 Conclusions 

It could be shown that the Social Force Model as proposed in 1995 (elliptical specification I) and 2000 (circular 
specification) in a special and one-dimensional case and around equilibrium reduces to the equations of the 
damped harmonic oscillator. This implies that indeed one has to be careful not to choose parameters with which 
the pedestrians’ trajectories yield unrealistic results. However, at the same time it means that there are parts in 
parameter space in which there are not only just small oscillations, but where there are exactly no oscillations. 
A look at parameter values found in literature for the circular specification shows that parameters deemed to 
yield realistic results in certain calibration scenarios do or may produce oscillating behavior unless the desired 
walking speed(s) are set to rather small values. 

The equations of the Social Force Model as of 2007 (elliptical specification II) are not as easily treated 
analytically. Still in a discussion of its equations it could be argued that - compared to the circular specification 
- oscillations are clearly suppressed. Elliptical specification II from 2007 therefore appears to be superior to the 
two preceding variants also for the reasons discussed in this paper (this was found already in the paper from 2007 
but for different reasons). It is therefore a good idea to either simply use elliptical specification II or combine it 
with one of the earlier variants (e.g. by simply adding the forces) to reduce the risk of oscillations. 

The phenomenon of oscillations was used to verify a specific computer implementation of the Social Force 
Model. It was possible to show that it reproduces the expected results. This comparison can be seen as an 
attempt to falsify either of the two - theoretical expectations and software implementation - and the attempt 
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failed, no potential issues could be found. The method can be applied to verify any implementation of the 
Social Force Model variants of 1995 or 2000 and it is generally an example of model specific verification. The 
tests carried out in this contribution are just examples. The phenomenon of oscillations bears the potential to 
formulate further tests. 
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A Proof that a pedestrian heading for a destination point never 
conies to rest 


In this section we show that 

n 

i=l 

diverges for n ^ oo. 

Eirst step: We proof that an+i < Oin- we begin with ao = 2 and with equation (24) 


Eor all 1 < z < 2 it holds that 


and for all—l/e< 2:<0 


Therefore for all 1 < z < 2 it holds that 


Q^n+l — 2 + Wo ( —) 


- < 0 


Wo{z) < 0 

§-W„iz) > 0 


Wo{-ze-^) < 0 

f^Woi-ze-^) > 0 

Erom that and with ao = 2 it follows that 

Oil — 2 Wo^ — OLqC ‘^°) < OLQ 


(67) 


( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 


Erom equation (24) it follows that 

a, - an+i = Wo(-an-ie-"-i) - Wo(-ane-"-) (76) 

and from it that if an-i > c^n also an > c^n+i- See figure 8 for a plot of 2 + Wo and the visualization 

of the evolution of the an which allows to see this first step very easily. 

Second step: We verify that Atn > 0 (physically this is obvious, mathematically it needs to be shown; this 
step in effect is a check on potential errors done before): Since 

an = -Wi (-ane-“'‘) (77) 

and because for — 1/e < z < 0 it holds that 

Wo(^)>Wi(z) (78) 
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Figure 8: Plot of 2 + Wq {—xe for 1 < a < 2 and a visualization of the evolution of the basically created 
with [20]. 


it results for the Atn according to equation (25) 

=an + Wo (-ane““”) = Wo (-aue““’*) - Wi > 0 (79) 

Third step: we proof that Atn+i < Atn- Equations (25) and (24) read 

= + Wo(-ane““'‘) (80) 

= 2 +Wo(-au-ie““’-^). (81) 

From them it follows that 

= a,+i + an - 2 (82) 

T 

thus the Atn show the same behavior as the an and with On+i < On also At^+i < Atn- 

Fourth step: if Atn+i/Atn for large n would approach a value smaller than 1 the sum would converge, only 
if it approaches 1 it may diverge. From (82) it follows that 


Atn-\-l 

T 

Oi-n 


Atn-\-l ^ an — 1 OLn-\-l 

Atn O^n — l 2 

To compute the limit for this for large n we have to expand Wo{—ane~°^'^) into a series. 
For convenience we write 

f3n — OLn 1 

knowing that this is ^n — UnjvQ and that these /3n will approach 0 with increasing n. 

By computing the right side limit towards 0 of 

Wo(-a„e-“") = Wo(-(l + 


(83) 


(84) 


(85) 


and its derivatives one can write 


Wo(-(l + /?n)e“(^+'’">) 


— 1 + /?n — g/?n + g/^n “ 405^" C>(/?n) 


compare figure 9. 

Considering only terms to second order 


Pn 


/3n-l 

/3n-l 


2 

3 

4 
3 


fn-l 

fn -1 


Using these approximations in equation (83) gives for large n 


Atn-\-l 

Atn 

which approaches 1 with n ^ oo where /3n-i ^ 0. 


2/3n-l 

3 - ^n-l 


( 86 ) 


(87) 

( 88 ) 


( 89 ) 
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Thus, like with the harmonic series hm the ratio of subsequent terms approaches 1. We note that equation 
(82) implies —= j3n+i + /3n and therefore 

oo 

=ao-l + 2Y,(3n (90) 

n = l 

the sum of the Atn depends trivially on the sum of the fin and especially their convergence behavior is the same. 
In the last step we prove that if fin > hm then also fin-\-i > hm+i and therefore beyond some n, m the harmonic 
series is a lower estimate of the fin implying that with the harmonic series also the sum of fin diverges. 

Assume that for some n, fin can be written as 


Pn=- (91) 

m 

with m G N, m >> 1 and I < q < ml{m — 1), i.e. fin is placed between two elements of the harmonic series. 
This requires just that there is some 0 < fin < I- This is obviously the case. 

So, can there be a fin+i < l/(m + 1)? 


fin-\-l 


1 

m + 1 


Q_ _ 2 _^ _ 1 

m 3 m + 1 


(92) 


A lower estimate of the right side is if at the positive term q is replaced by its minimum value q = 1 and at the 
negative term by its upper limit q — m/m—\\ 


fin+l 


1 

m + 1 


12 1 1 
m 3 (m — 1)2 m + 1 
mfi — Sm + 3 
^ 3m(?Ti — 1)2(?ti + 1) 

> 0V(m > 8) 


(93) 

(94) 

(95) 


Thus, if there is a fino < 1/8 - it is, as can be seen in table 1) - it will hold for all fin>no that if fin > hm then 
also fin+i > hm+i- Therefore the harmonic series is a lower estimate for the series of fin- Thus the series of the 
fin diverges and with it the series of the Atn- 



Figure 9: Plot comparing W{—{l^y) exp( —(1 + ^))) and its approximation [21]. 
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